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We report on Monte Carlo simulations focused on elucidating the phase structure of a SU (2) 
gauge theory containing Nt Dirac fermion flavors transforming in the fundamental representation 
of the group and interacting through an additional chirally invariant four fermion term. Pairs 
of physical flavors are implemented using the two tastes present in a reduced staggered fermion 
formulation of the theory with the Yukawa interactions necessary for generating the four fermion 
term preserving the usual shift symmetries. We observe a crossover in the behavior of the chiral 
condensate for strong four fermi coupling associated with the generation of a dynamical mass 
for the fermions. At weak gauge coupling this crossover is consistent with the usual continuous 
phase transition seen in the pure (ungauged) NJL model. However, if the gauge coupling is strong 
enough to cause confinement we observe a much more rapid crossover in the chiral condensate 
consistent with a first order phase transition 



The 30 International Symposium on Lattice Field Theory - Lattice 2012, 
June 24-29, 2012 
Cairns, Australia 



* Speaker. 

"'"Work supported in part by DOE grant DE-FG02-85ER40237 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Four fermion operators 



Simon Catterall 



1. Introduction 

Elucidating the nature of the electroweak symmetry breaking sector of the Standard Model 
(SM) is the main goal of the Large Hadron Collider currently running at CERN. It is widely be- 
lieved that the simplest scenario involving a single scalar Higgs field is untenable due to the fine 
tuning and triviality problems which arise in scalar field theories. One natural solution to these 
problems can be found by assuming that the Higgs sector in the Standard Model arises as an 
effective field theory describing the dynamics of a composite field arising from strongly bound 
fermion-antifermion pairs. These models are generically called technicolor theories. 

However, to obtain fermion masses in these scenarios requires additional model building, as 
in extended technicolor models ^, ||, fij] and models of top-condensation [||, |6[ f^, ||]. In the latter 
models four-fermion interactions drive the formation and condensation of a scalar top-anti-top 
bound state which plays the role of the Higgs at low energies. 

Our motivation in this paper is to study how the inclusion of such four fermion interactions 
may influence the phase structure and low energy behavior of non-abelian gauge theories in general. 
Specifically we have examined a model with both gauge interactions and a chirally invariant four 
fermi interaction - a model known in the literature as the gauged NJL model [^]. 

The focus of the current work is to explore the phase diagram when fermions are charged 
under a non-abelian gauge group. Indeed, arguments have been given in the continuum that the 
gauged NJL model may exhibit different critical behavior at the boundary between the symmetric 
and broken phases 1 corresponding to the appearance of a line of new fixed points associated with a 
mass anomalous dimension varying in the range 1 < < 2 || [TOQ . The evidence for this behavior 
derives from calculations utilizing the ladder approximation in Landau gauge to the Schwinger- 
Dyson equations. A primary goal of the current study was to use lattice simulation to check the 
validity of these conclusions and specifically to search for qualitatively new critical behavior in the 
gauged model as compared to the pure NJL theory. While we will present results that indicate that 
the phase structure of the gauged NJL model is indeed different from pure NJL, we shall argue that 
our results are not consistent with the presence of any new fixed points in the theory. 

In the work reported here and described in detail in [jll|] we have concentrated on the four flavor 
theory corresponding to two copies of the basic Dirac doublet used in the lattice construction. The 
four flavor theory is expected to be chirally broken and confining at zero four fermi coupling and 
is free from sign problems for gauge group SU (2). Understanding the effects of the four fermion 
term in this theory can then serve as a benchmark for future studies of theories which, for zero 
four fermi coupling, lie near or inside the conformal window. In the latter case the addition of 
a four fermion term will break conformal invariance but in principle that breaking may be made 
arbitrarily small by tuning the four fermi coupling. It is entirely possible that the phase diagrams of 
such conformal or walking theories in the presence of four fermi terms may exhibit very different 
features than those seen for a confining gauge theory. 



Notice that the appearance of a true phase transition in the gauged NJL models depends on the approximation that 
we can neglect the running of the gauge coupling 
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2. Details of the model 

We will consider a model which consists of Nf/2 doublets of gauged massless Dirac fermions 
in the fundamental representation of an SU (2) gauge group and incorporating an SU (2)l X SU (2)r 
chirally invariant four fermi interaction. The action for a single doublet takes the form 



d 4 xy(i$ -4)V^- ^[(W) 2 + (W5T» 



1 

2g 



2 7>[F MV FH, (2.1) 



where G is the four-fermi coupling, g the usual gauge coupling and r a , a = 1 ... 3 are the generators 
of the SU (2) flavour group. 

This action may be discretized using the (reduced) staggered fermion formalism with the result 

S = £ X T {x) %{x) fo M (*) + G0 M (x)e(x)^(x)]. (2.2) 

where rj^(x), ^(x) and e(x) are the usual staggered fermion phases, <f>(x) = j^Y.iA{ x ~ h) the 
average of the scalar field over the hypercube [12, 13] and the gauge field acting on the reduced 
staggered fermions takes the form: 

%(*) = ^[1 + e(*)] U^x) + ~[1 -e(x)} u;(x). (2.3) 

Clearly the theory is invariant under the U(l) symmetry %(x) — > e' ae ( x 'X(x) which is to be 
interpreted as the U(l) symmetry corresponding to fermion number. More interestingly it is also 
invariant under certain shift symmetries given by 

Xix) -> $ p (x)x(x + p), (2.4) 
Upix) ^U;(x + p), (2.5) 
<^(x) -► (-l) 5 ^(x + p). (2.6) 

These shift symmetries correspond to a discrete subgroup of the continuum axial flavor transfor- 
mations which act on the matrix field *P according to 

¥ -> js^Yp (2-7) 
Notice that no single site mass term is allowed in this model. 



3. Numerical results 

We have used the RHMC algorithm to simulate the lattice theory with a standard Wilson gauge 
action being employed for the gauge fields. Upon integration over the basic fermion doublet we 
obtain a Pfaffian Pf (M(U)) depending on the gauge field 2 . The required pseudofermion weight for 
Nf flavors is then Pf(M) N ^ 2 . The pseudoreal character of SU (2) allows us to show that the Pfaffian 

2 Note that the fermion operator appearing in eqn. is antisymmetric 



3 



Four fermion operators 



Simon Catterall 



is purely real and so we are guaranteed to have no sign problem if we use multiples of four flavors 
corresponding to a pseudofermion operator of the form (M^M) s - . The results in this paper are 
devoted to the case Nf = 4. We have utilized a variety of lattice sizes: 4 4 , 6 4 , 8 4 and 8 3 x 16 and 
a range of gauge couplings 1.8 < j6 = 4/g 2 < 10.0. To determine where the pure gauge theory is 




Figure 1: Polyakov loop vs j3 at G = 0.1 for four flavours 



strongly coupled and confining we have examined the average Polyakov line as /3 varies holding the 
four fermi coupling fixed at G = 0.1. This is shown in figure [[]. We see a strong crossover between 
a confining regime for small jS to a deconfined regime at large /3 . The crossover coupling is volume 
dependent and takes the value of j3 c ~ 2.4 for lattices of size L = 8. For j8 < 1 .8 the plaquette drops 
below 0.5 which we take as indicative of the presence of strong lattice spacing artifacts and so we 
have confined our simulations to larger values of j3. We have set the fermion mass to zero in all of 
our work so that our lattice action possesses the series of exact chiral symmetries discussed earlier. 

One of the primary observables used in this study is the chiral condensate which is computed 
from the gauge invariant one link mass operator 

X (x) ( % (x)x (x + ep) + &Ji(x-ep)x{x-ep))e (x) ^ (x) (3.1) 

In Figure || we show a plot of the absolute value of the condensate at a variety of gauge 
couplings j8 on 8 4 lattices. Notice the rather smooth transition between symmetric and broken 
phases around G ~ 0.9 for /3 = 10. This is consistent with earlier work using sixteen flavors of 
naive fermion reported in jl4] ] which identified a line of second order phase transitions in this 
region of parameter space. It also agrees with the behavior seen in previous simulations using 



conventional staggered quarks [15] 



This behavior should be contrasted with the behavior of the condensate for strong gauge cou- 
pling j8 < 2.4. Here a very sharp transition can be seen reminiscent of a first order phase transition. 
In Figure |3] we highlight this by showing a plot of the condensate versus four fermi coupling at 
the single gauge coupling j8 = 2.0 for a range of different lattice sizes. The chiral condensate is 
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Figure 2: (%%) vs G for varying /3 for the 8 4 lattice with Nf = 4. 
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Figure 3: (%X) vs G at j3 = 2.0 for lattices 4 4 , 6 4 and 8 4 with N f = 4. 



now non-zero even for small four fermi coupling and shows no strong dependence on the volume 
consistent with spontaneous chiral symmetry breaking in the pure gauge theory. However, it jumps 
abruptly to much larger values when the four fermi coupling exceeds some critical value. This 
crossover or transition is markedly discontinuous in character - reminiscent of a first order phase 
transition. Indeed, while the position of the phase transition is only weakly volume dependent it 
appears to get sharper with increasing volume. 

What seems clear is that the second order transition seen in the pure NJL model is no longer 
present when the gauge coupling is strong. In the next section we will argue that this is to be ex- 
pected - in the gauged model one can no longer send the fermion mass to zero by adjusting the 
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four fermi coupling since it receives a contribution from gauge mediated chiral symmetry break- 
ing. Indeed the measured one link chiral condensate operator is not an order parameter for such 
a transition since we observe it to be non-zero for all G. Notice however that we see no sign that 
this condensate depends on the gauge coupling j3 in the confining regime at small G. This is qual- 
itatively different from the behavior of regular staggered quarks and we attribute it to the fact that 
the reduced formalism does not allow for a single site mass term or an exact continuous chiral 
symmetry. Thus the spontaneous breaking of the residual discrete lattice chiral symmetry by gauge 
interactions will not be signaled by a light Goldstone pion and the measured condensate will re- 
ceive contributions only from massive states. The transition we observe is probably best thought of 
as a crossover phenomenon corresponding to the sudden onset of a new mechanism for dynamical 
mass generation due to the strong four fermi interactions. 



4. Summary 

In this paper we have conducted numerical simulations of the gauged NJL model for four 
flavors of Dirac fermion in the fundamental representation of the SU (2) gauge group. We have 
employed a reduced staggered fermion discretization scheme which allows us to maintain an exact 
subgroup of the continuum chiral symmetries. 

We have examined the model for a variety of values for lattices size, gauge coupling, and 
four fermi interaction strength. In the NJL limit /3 — >■ °° we find evidence for a continuous phase 
transition for G ~ 1 corresponding to the expected spontaneous breaking of chiral symmetry. How- 
ever, for gauge couplings that generate a non-zero chiral condensate even for G = this transition 
or crossover appears much sharper and there is no evidence of critical fluctuations in the chiral 
condensate. 

Thus our results are consistent with the idea that the second order phase transition which exists 
in the pure NJL theory (/3 = °°) survives at weak gauge coupling. However our results indicate that 
any continuous transition ends if the gauge coupling becomes strong enough to cause confinement. 
In this case we do however see evidence of additional dynamical mass generation for sufficiently 
large four fermi coupling associated with an observed rapid crossover in the chiral condensate and 
a possible first order phase transition. 

The fact that we find the condensate non-zero and constant for strong gauge coupling and 
G < G c shows that the chiral symmetry of the theory is already broken as expected for SU (2) with 
Nf = 4 flavors. This breaking of chiral symmetry due to the gauge interactions is accompanied by 
the generation of a non-zero fermion mass even for small four fermi coupling. Notice that this type 
of scenario is actually true of top quark condensate models in which the strong QCD interactions 
are already expected to break chiral symmetry independent of a four fermion top quark operator. 
The magnitude of this residual fermion mass is not controlled by the four fermi coupling and cannot 
to sent to zero by tuning the four fermi coupling - there can be no continuous phase transition in the 
system as we increase the four fermi coupling - rather the condensate becomes strongly enhanced 
for large G. 
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